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Abstract 

We examine the weak noise limit of an overdamped dissipative system within 
a semiclassical description and show how quantization influences the growth 
and decay of fluctuations of the thermally equilibrated systems. We trace 
its origin in a semiclassical counterpart of the generalized potential for the 
dissipative system. 
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I. INTRODUCTION 



One of the important issues in nonequilibrium phenomena in the macroscopic nonlinear 
systems is to understand the interplay of nonlinearity of the system and the fluctuations 
of its environment. The problem is fairly general in the context of chemical reactions [|l|, 
breakdown of electronic devices [EJ, phase transitions 0] etc. The essential description of 
the physical situation rests on the Fokker-Planck equations for the probability distribution 
functions of the relevant variables of the dynamics. In the weak noise limit the fluctuations 
have been described |3|-|| by appropriate auxiliary Hamiltonian or path integral methods. 
The theoretical results have been corroborated by remarkable experiments on fluctuations 
using analogue electronic circuits || , which allow the phase space trajectories of fluctuations 
to be observed directly in a precise manner. These studies have enriched our understanding 
in several theoretical issues, e. g. , the symmetry between the growth and the decay of 
classical fluctuations in equilibrium and its breakdown under nonequilibrium conditions |3j 
, the existence of a nonequilibrium potential of a dissipative system |5]|| etc. 

It is the purpose of this paper to extend the theory to the semi classical context. The 
quantization of the system itself adds a new dimension to the interplay of nonlinearity and 
stochasticity in a dissipative system. To make a fair comparison with classical theory we 
adopt the Wigner's phase space distribution function of c-number variables. The weak noise 
limit can then be appropriately employed to develop an auxiliary Hamiltonian formulation 
at the semiclassical level in terms of these phase space variables. This allows us to realize 
the existence of an optimal force of purely quantum origin derivable from the fluctuating 
field and relate it to the momentum of the auxiliary Hamiltonian. The quantum correction 
also makes its presence felt in the growth and decay of fluctuations of thermally equilibrated 
semiclassical systems keeping the symmetry preserved. 

The outline of the paper is as follows: In Sec. II we introduce the general aspects of 
dynamics of dissipative quantum system in terms of the Wigner phase space function. In 
Sec. Ill we consider the weak noise and semiclassical limit under overdamped condition and 
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take resort to the well-known auxiliary Hamiltonian description. The quantum part of the 
optimal force is then explicitly derived. The symmetry between the growth and decay of 
fluctuations in a thermalized quantum system is discussed in Sec. IV. The existence of a 
semiclassical contribution to the potential in a dissipative system is then shown in Sec. V. 
The paper is concluded in Sec. VI. 



II. QUANTUM DYNAMICS OF A DISSIPATIVE SYSTEM 

We consider a dynamical system characterized by a potential V(x) coupled to an envi- 
ronment. Evolution of such an open quantum system has been studied over the last several 
decades under a variety of reasonable assumptions [[5]-§| . Specifically interesting is the semi- 
classical limit of an Ohmic environment. The dissipative time evolution of the Wigner 
distribution function W(x,p, t) for the system with unit mass (m = 1) under the potential 
V(x) can be described by []l0H l~2"ll 



dW _ dW dVdW v h 2n (-l) n d 2n+1 V d 2n+l W dpW d 2 W 

where 7 and D are the dissipation constant and the diffusion coefficient, respectively, x and p 
are c-number co-ordinate and momentum variables. The drift term is a direct consequence || 
of the existence of 7-dependent term in the imaginary part of the exponent in the expression 
for the propagator for the density operator in the Feynman-Veron theory and has been shown 
to be responsible for appearance of a damping force in the classical equation of motion for 
the Brownian particle to ensure quantum-classical correspondence. 7 and D are related by 
the fluctuation-dissipation relation, D = ^Uuj coth (in the semiclassical limit D = jkT 
). u is the renormalised linear frequency of the nonlinear system. The quantum correction to 
classical Liouville motion is contained in the ^-containing terms in the sum. e is a parameter 
(whose value is 1) which is kept in the equation for bookkeeping the Wigner correction term 
in our further analysis. We put e = 1 at the end of calculation. 

Eq.(l) had been used earlier in several occasions. For example, Zurek and Paz [IIJ and 



others [|l2j have studied some interesting aspects of quantum-classical correspondence in 



relation to decoherence and chaos. Based on this equation and its variant chaotic dissipative 



systems has been studied. |IT]fl3| -|H)|| . The equation also yields the simplest leading order 
quantum correction term to classical Kramers' rate |fTE| . The primary reason for choosing 
Eq.(3) as our starting point is that it reaches the correct classical limit when h — > so 
that D becomes a thermal diffusion coefficient in the high temperature limit and the 

Wigner function reduces to the corresponding classical phase space distribution function and 
we recover the Kramers' equation which describes classical Brownian motion of a particle in 
phase space. 

III. WEAK NOISE AND SEMICLASSICAL LIMIT OF QUANTUM DISSIPATIVE 

DYNAMICS: 

The weak noise limit of a dissipative system within a semiclassical description can be 
conveniently described by a "WKB-like" ansatz (we refer to "WKB-like" since we are con- 
sidering more than one dimension. Traditionally WKB refer to one dimension only) of the 
Eq.(l) for the Wigner function of the form 

W{x,p,t) = Z{x,t)exp{-^)exp{-£-) . (2) 

where D\ — The weak noise limit is defined || as D\ — ► and semiclassical limit 
refers to (K —* 0). Z(x, t) is a prefactor and s(x,p, t) is a classical action which is a function 
of c-number variables x and p , satisfying the following Hamilton- Jacobi equation, 

ds ^ ds y>^ s ^ 2^ s y 

dt dx dp dp dp 

4l 2 2n {2n)\ dx 2n+1 dp V } 

The derivation of Eq. (3) is based on the following consideration. Since in the weak noise 
limit Di is the relevant small parameter one obtains with ansatz (2) in leading order a term 
proportional to h n yj^^j which is not balanced by any other term of the same order 
Df (2n+1) . This is because the highest derivative in Eq.(l) does not have a factor scaling 



with the corresponding power of D\. The successive terms next to the leading order are also 
/i-containing terms. All of these terms vanishes in the semiclassical limit (% — > 0). Therefore 
the leading order term that remains gives rise to Eq.(3). It is thus obvious that an ansatz 
(2) with h finite is not feasible. The semiclassical limit h — > is a necessary requirement for 
the validity of ansatz (2). The above equation can be solved by integrating the Hamiltonian 
equations of motion, 

x = p 

X = P-jX 



p 



-l) 3n+1 d ( ^d 2n+1 V 



x 2n 



^ 2 2n (2n)\ dx \ dx 2n+1 
which are derived from the following effective Hamiltonian H e ff 



X (4) 



H eff =pP- V'X - jXp + u 2 X 2 



2 2 "(2n)! dx 2n+lX ' ^ 

Here we have put ^ = P and |^ = X. The introduction of additional degree-of-freedom 
by incorporating the auxiliary momentum (P) and co-ordinate (X) makes the system an 
effectively two-degree-of-freedom system. The origin of these two variables is the thermal 
fluctuations of the environment ||. The auxiliary Hamiltonian H e ff is not to be confused 
with the microscopic Hamiltonian comprising the system, the bath and their coupling. Thus 
the phase space trajectories concern fluctuations of the c-number variables. 

Under overdamped condition (it <C 7x5 X <C jX) Eqs. (4) can be easily reduced to the 
following form ; 

2u 2 X 



X = K(x) 
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where 



K(x) = - 

7 



v (-1) 3 "+V" d 2n+1 V 



(7) 



It is easy to recognize the quantity 2uJ ^ x as a momentum and redefine it as p r . Therefore 
Eqs(6) may be rewritten as 

X = K(x) + p r 

dK(x) 

Pr = ~^Pr ■ (8) 

So the overdamped motion is described by the following effective Hamiltonian, 

2 

Hod = j + K{x)p T . (9) 

The above auxiliary Hamiltonian description (8, 9) is isomorphic in form to that of 
Luchinsky and McClintock |§, who had considered an overdamped classical Brownian mo- 
tion in a force field K(x), driven by a weak white noise ((t) whose intensity D\ <C 1 as 

x = K(x) + ((t), (C) = , <C(*)C(0)> = D 1 S(t) . (10) 
Equivalently the corresponding Fokker-Planck equation for the probability density P(x, t) 

is 

dP(x, t) 8K(x) D l d 2 P{x,t) 

dt ~ dx ^ x ^ + 2 dx 2 ' 1 J 

The large fluctuations of scale 3> \fD% can therefore be treated in the weak noise limit 
Di — > by "WKB like" approximation of the Fokker-Planck equation (11) in the 

form 

P(x,f) = y(x,f)expP^)] . (12) 

Here y(x, t) is the prefactor and </>(#, t) is a "classical" action describing a Hamiltonian- 
Jacobi equation which can be solved by solving Hamilton' s equation (8) with p r = |^ as 
the momentum for auxiliary system. 
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The important distinguishing feature of the above description in which the system is 
treated semiclassically is the structure of K(x) which is given by equation(7) and comprises 
of two terms ; 

K{x) = K d + K semi {x) (13) 

where 

K* - _m (14) 

7 

is derivable from purely classical potential V(x) and K semi (x) does not explicitly involve 

% 

K semi (x) = - E 2 2 "(2n)! dx 2 ^ 1 (15) 
originates from the nonlinearity of the potential V(x) and quantum nature of the system. 
The quantum contribution to K(x) is therefore likely to influence both the fluctuational and 
the relaxational paths of the dynamics and is also responsible for the existence of a potential. 
Our objective is now to explore these aspects in the following two sections. 



IV. LARGE FLUCTUATIONS IN EQUILIBRATED SEMICLASSICAL SYSTEMS 

In the thermally equilibrated systems a typical large fluctuation of the variable x implies 
a temporary departure from its stable state, x s to some remote state Xf. This is followed 
by a return to x s as a result of relaxation in the absence of fluctuations p r . A nonzero 
value of p r which results from fluctuations due to surrounding drives the system away from 
x s along a set of trajectories which form the unstable invariant manifold and define the 
so called fluctuational paths. On the other hand the system relaxes along the relaxational 
return path to x s under the condition p r = 0, which form stable invariant manifold. The 
latter condition implies x = K(x). In each case the trajectory represents the optimal paths 
along which the system is expected to move with overwhelming probability. Luchinsky and 
McClintock || have studied these paths in analog electronic circuits and demonstrated the 
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growth and the decay of classical fluctuations in equilibrium. We extend this analysis to 
semiclassical domain using the same model potential, 

V(x) = \x' - \x" . (16) 

The quantum contribution to the growth and the decay of fluctuations can be under- 
stood by recognizing the K semi (x) term in the dynamics (8). In Fig.l we compare both 
the deterministic fluctuational and relaxational (optimal) paths for quantum and classical 
thermally equilibrated systems. It is important to note that the maximum possible ampli- 
tude of large fluctuations is almost double for the quantum system compared to that for the 
corresponding the classical system. This is due to the addition of the nonlinear force term 
of quantum origin, K semi (x) in the c-number equation (8), which is shown to be derivable 
from the fluctuating field, and is related to the momentum of the auxiliary Hamiltonian. 

Before leaving this section we make a brief remark on the thermally nonequilibrated 
systems. Since the detailed balance is not operative here, the optimal path to a given state 
not just the time- reversed dynamical path along which the system moves from this state 
to the stable state in absence of fluctuations p r . Thus for the driven system, for example, 
the pattern of optimal path is generically different from that for the thermally equilibrated 
systems. It may display singularities whose topological manifestations as caustics, switching 
line, cusps etc have been thoroughly studied for classical systems. We believe that K semi (x, t) 
where t signifies the driving by a periodic force in V(x, t) is likely to play an important role 
in their quantum counterparts. 

V. EXISTENCE OF A POTENTIAL FOR DISSIPATIVE SEMICLASSICAL 

SYSTEM 

In a significant analysis Graham and coworkers had examined the condition for 
existence of potential for classical dissipative systems. We now extend this analysis to the 
present semiclassical context. 
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The general criterion for a dissipative dynamical system described by autonomous equa- 
tions of the form 



x v = K v (x) (17) 

to have a potential (f>(x) with respect to Q v ^ (positive, semidefinite symmetric matrix, 
which are considered to be the matrix of transport coefficients) if there exists a single- valued 
continuously differentiable and globally defined function (f>(x), bounded from below which is 
stationary in the limit sets of Eq.(17) and which satisfies 

IC(x) = -\cr»^ + r*(x) (18) 

with 

r^ = 0. (19) 

dx v v ; 

Here for simplicity Q Vfx is assumed to be independent of x. The first and the second 
terms of Eq.(18) correspond to irreversible and reversible part, respectively. 

The stochastic process x(D 1 ,t) which involve Eq.(17) and a symmetric non-negative 
matrix Q v ^ is governed by the Fokker-Planck equation for probability distribution function 

P(x,t) 

For Di — the above description reduces to (17). For D\ ^ the steady state distribu- 
tion defines the function <p(x, t) by 

P(x, D u t -> oo) = N(D 1 ) exp[-^^] (21) 

N is the normalization constant. If <f>{x) = lim^^o 4>( x i Di) is a single-valued, continu- 
ously differentiable function bounded from below it satisfies 
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Eq. (22) is equivalent to Eqs. (18, 19). Interpreting Eq.(22) as usual as a Hamilton- 
Jacobi equation by defining (p(x) as action and J^- = P v as the momentum conjugate to x u , 
one can construct the auxiliary Hamiltonian 

H(x,p) = ^Q^VvVv + K v (x) Vu . (23) 

Graham and co-workers |5]|| have argued that a potential can exist with equation (17) if 
the above Hamiltonian is integrable for H = 0, because the condition implies that there exist 
a smooth separatrix which connect smoothly the stable and unstable manifolds emanating 
from the hyperbolic fixed points of the dynamical system. 

We now turn back to our dissipative semiclassical system described by Eqs(8) and (9) 
where K(x) is defined by Eq.(7). Recognizing Eq.(9) as Eq.(23) for a one-degree-of- freedom 
system we identify 

Q = l 

K(x) = K d + K semi . (24) 
The potential function <j)(x) can therefore be calculated from Eq.(9) with H = as (since 

_ d<t>(x) \ 
Vr dx I 

<f)(x) = -2 J K(x)dx . (25) 

The above expression can be made more explicit if we make use of Eq.(7) in (25). We 
obtain 

<f>(x) = (f)d(x) + <p S emi{x) (26) 



where 



2 r _ 2 



V\x)dx = -V{x) (27) 



7 J 7 



and 



2 _ r x 2n (-l) 3n+1 d 2n+l V , 
( x ) = Y ; - '— - - - , 1 dx 28 



10 



The existence of a potential for a dissipative, semiclassical dynamical system is thus 
ascertained. The method essentially relies on a dynamical definition p r as a derivative of 
the potential (f>{x) in a system described by an overdamped quantum Markov process in the 
weak noise and semiclassical limit. As elaborated earlier in Sec. Ill p r has a statistical origin 
which drives the system away from its stable state x s to a preassigned remote state Xf from 
which the system relaxes in absence of p r . It is thus important to realize that both the 
dynamical and statistical notions are kept intact in the quantum treatment. 

VI. CONCLUSIONS 

Keeping in view of the quantum nature of the system in terms of the Wigner's phase 
space function we examine the semiclassical dynamics of a dissipative system in an Ohmic 
environment. The weak noise limit of the stochastic process then allows us to capture 
the essential features of the dynamics within the framework of an auxiliary Hamiltonian 
description at the semiclassical level. Our results are summarized as follows : 

(i) The Wigner's quantum correction to classical Liouville equation gives rise to an optimal 
force in addition to usual the classical force term. This quantum optimal force is essentially 
a result of an interplay of nonlinearity of the system and the thermal fluctuations of its 
environment and is derivable in terms of an auxiliary Hamiltonian description. 

(ii) This term is also responsible for modification of growth and decay of large fluctuations 
from equilibrium for the appropriately thermalized quantum system (compared to its clas- 
sical counterpart). The symmetry of the fluctuational and the relaxational paths, signifying 
the detailed balance, however, as expected is kept preserved. 

(iii) The quantum correction term implies the existence of a potential for the dissipative 
semiclassical system. 
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Since the fluctuational and the relaxational paths have been experimentally demonstrated 
as a part of physical reality by analogue experiments || in the realm of large fluctuations, 
we believe that the essential modification of the integrable, nonintegrable and the singular 
topological features of the dynamics due to semiclassical correction might be relevant in 
several contexts. We hope to address some of these issues in a future communication. 
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FIGURES 

FIG. 1. A plot of system co-ordinate x vs time t signifying the fmctuational and relaxational 
paths according to Eq.(8) for the model system with V(x) = ^x 4 — ^x 2 for (a) semiclassical case 
(e = 1) and (6) classical case (e = 0), (units arbitrary). 



14 



B. C. Bag and D. S. Ray 




